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Abstract 

A q-difference analogue of the Painleve III equation is considered. Its derivations, 
affine Wcyl group symmetry, and two kinds of special function type solutions are 
discussed. 



■ 1 Introduction 

G\ ' 

Importance of the Painleve equations is widely recognized in mathematics and mathemat- 
ical physics, since the solutions plays a role of nonlinear version of special functions. So 
much studies for the Painleve equations have been done from various points of view, such 
as classical solutions, asymptotics, geometric or algebraic structures and so on ^3 121- On 
the other hand, discrete analogue of the Painleve equations have been studied extensively 
after discovery of the singularity confinement property [4, which is a discrete analogue 
of the Painleve property. What matters in studies of the discrete Painleve equations was 
that they were in some sense too rich and looked almost uncontrollable. For example, 
number of known discrete Painleve equations is several ten's, while that for the Painleve 
equations is only six. 

The situation has been changed after Sakai's theory {22 . Sakai constructed "space of 
initial conditions", or defining manifold of a discrete Painleve equation called q-Pyi by 
blowing up P 2 at nine points. Based on this fact, Sakai gave a classification of discrete 
Painleve equations with the aid of theory of rational surfaces, and discussed the symme- 
tries. The Painleve equations appear naturally from the degeneration of defining manifold 
which is also interpreted as continuous limit. This theory gives us clear view to both 
discrete and continuous Painleve equations. It also suggests the list of discrete Painleve 
equations to be studied further. 
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In this paper, we shall work with the following (/-difference equation 



fo 



q 2N + l c 2 1 + aoq n fo 

/o/i a q n + f ' 
q 2N + 1 c 2 a iq - n+v + h 
Mi 1 + aig-"+"/i 



(1.1) 



where fa = fi(n;v,N) (i = 0,1) are dependent variables, n £ Z is the independent 
variable, v,N € Z are parameters, and (/, do, ai, c are constants. Moreover, /j and /, 
denote fi(n + l;u, N) and fi(n — 1; z/, iV), respectively. We call equation (jl.lj) a (/-difference 
Painleve III equation(g-Pm). Since (/-Pm appears in "Mul.6" in Sakai's classification, it 
is expected that g-Pm is a "good" equation. 

The purpose of our work is to establish and to give explicit description for the following 
properties of g-Pm (jl.ljl : 

• (/-Pm admits a continuous limit to the Painleve III equation (Phi)) 



where a, (3 are parameters. 

• (/-Pm is derived from the discrete-time Relativistic Toda lattice equation [23] , while 
Pni is derived from the relativistic Toda lattice equation |21| . 

• (/-Pm can be regarded as a discrete dynamical system on the root lattice of type 
A2 x Ai j221 15]. As a consequence, it admits symmetry of (extended) affine Weyl 
group of type x as the group of Backlund transformations. 

• q-Pin admits two classes of classical solutions, namely, the solutions which are re- 
ducible to discrete Riccati equation, and rational solutions. There are two kinds of 
Riccati type solutions, one of which is expressed by the Jackson's (modified) g-Bessel 
functions, while the Riccati type solution of Pm is expressed by the (modified) Bessel 
functions [T7] . 

• The rational solutions of Pni (|l-2j) are expressed as ratio of some special polynomials 
(Umemura polynomials). They admit determinant formula of Jacobi-Trudi type 
whose entries are given by the Laguerre polynomials jlUj . The rational solutions of 
g-Pni are characterized by similar special polynomials which also admit determinant 
formula. 

In this paper, being the first half of our work, we discuss the derivation, symmetry, 
and particular solutions of Riccati type. Rational solutions will be discussed in the next 



This paper is organized as follows. In Section 2 we derive Pni and q-Pm from the rela- 
tivistic Toda lattice and discrete-time relativistic Toda lattice, respectively. Moreover, we 
also derive g-Pni from the birational representation of affine Weyl group of type 
which was presented in [5J. In Section 3, we construct two kinds of particular solutions of 
Riccati type, one of which is expressed by Jackson's modified (/-Bessel functions. We also 
construct determinant formula for these solutions. 




7 = -5 = 4, 



(1.2) 
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Remark 1. 1. Originally q-Pm was derived in as "asymmetric generalization" 
of g-difference Painleve II equation (q-Pu) [T5] . 

_ c x n + (A n . . 

1 + C, A x n 

where £, £ are parameters, and A is a constant. In fact, putting 
q = A -2 , i/ = 0, ai = — -, 

x 2 n = /i(n; 0, AT), x 2n+1 = / (n; 0, iV), (1.4) 
g-Pm CU) is reduced to g-Pn fTSj) with 

$ = A^+V, C = i (1-5) 

OoA 

2. The following ^-difference equation, 

ai(wn - ag n )(wn - &<? n ) c , 

Wn+1%-1 = 7 w ~7\ J I 1 - 6 ) 

(io n - c)(to w - d) 

is also identified as a ^-difference Painleve III equation |18M19| . Here, q is a constant and 
o, b, c, d are parameters. A class of particular solutions in terms of Jackson's g-Bessel 
functions is constructed in . It looks that equation Ijl.fijl has a different nature compared 
to our g-Pm (jl.l|) . and it may be natural to regard equation ()1.6|) as a degenerate case of 
the g-difference Painleve VI equation (g-Pyi) [71 122j. 

a 3 a 4 (z n+ i - 6ig n )(z n+ i - b 2 q n ) 
VnVn+l = ■ 



(Zn+i - h)(z n+ i - 64) 

_ bzb A (y n - aiq n )(y n - a 2 q n ) hb 2 _ a x a 2 , ■> 

z n z n +i — t T7 r , T~T ~ 1 ' v 1 -'] 

where a, and 6j (i = 1,2,3,4) are parameters. 

2 Derivations of q-P m 

2.1 Derivation of Pni from relativistic Toda lattice 

In this section we derive Pni (|1.2|) from the relativistic Toda lattice equation [211 123j , 

U n = d n (c n -c n ^ 
at 

-T x c n = c n (d n+ i + c n+ i - d n - Cn-i), n € Z. (2.1) 
dt 

We introduce the variables 14 and K n by 

KAT n V n -iK n -i 
Cn = ~-r; , d n = —K n _i -\ , (2.2) 
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which yields 

j t log V n = + V n -{) - (K n + V n ), 

—K n = K n V n - K n+1 V n+1 , n € Z. (2.3) 
We next impose two-periodicity on equation IJ2.MJI . Then we have, 

^-io g v = (K 1 + v 1 )-(K + v ), 

at 

^-logV 1 = (K + V )-(K 1 + V 1 ), 
at 

^-K = K V -K 1 V 1 , 
at 

^-K 1 = K 1 V 1 - K V . (2.4) 
at 

In order to derive Pm (|1.2|) . we introduce additional constants a n (n = 0, 1) in such a way 
that 

^ log V = {K x + Vi) - (K + Vq) + Q , 

^ log 14 = (K + Vb) - (X a + V X ) + ai, 
at 

-K Q = K V -K 1 V 1 , 
at 

-K 1 =K 1 V 1 -K V . (2.5) 

We note that integrability is still kept under this generalization in a sense that it admits 
the Lax pair. It is easy to see from equation 1)2. 5|) that we have Kq + K\ = const and 
4 log(VbVi) = ao + a±. We normalize them so that 

ao + ai = l, V V 1 =e t , K Q + K 1 =fo. (2.6) 

Under this normalization, one can write down the equations for Vq and K\. Eliminating 
Kx, we obtain 



d 2 



^ = W (^f) + V o ~ («o + Po)V 2 + (A) -ao + l)e* - ^, (2.7) 



dt 2 V \dt J ' " u ' ' ^ u " u ' ^ Vo 

or 

dVo 1 ^ \ 2 ldV V 2 /?o-ao + l 1 , . 

"d^^U^J -^ + ^ (y °- (ao + A))) + i V (2 - 8) 

where s = e , which is a version of Pm (named as "Pm'" in 17 ]). Equation (|2.8|) is 
rewritten into the ordinary form of Pm (|1.2[) with a = — 4(«o + A))j /3 = 4(— ao + A) + 1) ; 
7 = — (5 = 4 by putting s = x 2 and Vq = xv. 
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Remark 2. Without introducing additional constants a n , we obtain the equation for Vo 

as 

d?Vp _ 1 (dVpV 3 2 1 

Equation (|2.9j) is classified as "XII" in Gambier's classification and integrated in terms 
of the elliptic functions. Therefore, such generalization of the relativistic Toda lattice 
equation is crucial in order to derive Pin- Such situation can be seen for other Painleve 
equations, where they are derived by the similar generalization of proper soliton equa- 
tions p us nnnn. 

2.2 Derivation of q-Pm from discrete-time relativistic Toda lattice 

Let us next consider the discrete case. We start from the discrete-time relativistic Toda 
lattice equation 23 in the following form, 

dn 0"n— 1 d n — \ C n Q"n— 1 Cn— 1 

dn &n d n C n 

( - it 

a n = l + d n + ^±I, neZ, (2.10) 

where a n = a n (t), c n = c n (t), d n = d n {t) and denotes the value at t + 1. We note that 
for convenience we have inverted the space direction n of the original equation in 2Z^ ■ We 
introduce the variables K n = K n (t), V n = V n (t) by 

— = V n , a n = K n . (2.11) 

Then equation (|2,10|) is rewritten as 

Vn _ K n _i + V n -l 
~ K n + V n 

K n +1 _ 1 + Vn K n+ i + V n +l 
Kn ~ l + V n+1 Kn + V n 

We next impose two-periodicity on equation (|2.12|) . which yields 

Fp _ Ki + Vi V 1 _ K + Vp 

Vo Ko + Vo'_ Vi Ki + Vt' 

K l _ 1 + V K 1 + V l Kp _ 1 + Vi Kp + Vp 
Ko ~ 1+Fi Ko + Vo' K x ~ l+F Ki + V{ 



(2.12) 



(2.13) 



Similarly to the continuous case, we introduce additional constants a n (n = 0, 1) in such 
a way that 



Vo _ Ki + V 1 V 1 _ Kp + Vp 

Vo Ko + Vo^ _V X Ki + V^ 1 ^ ^ 

K x = 1 + Vp Ki + V 1 Kp = 1 + Vi Kp + Vp 

Kp ~ 1 + Vi Ko + Vp ' K x ~ l + F Ki + Vi ' 



(2.14) 
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We remark equation (|2.14| is still integrable in a sense that it passes the singularity 
confinement test J3]. It is easy to see from equation (j2.14j) that we have KqK\ = KqKi = 
const and Fo^i = aoaiVoVi. We normalize as 

a oai =q 2 , y y 1=c 2 g 2 *, K K 1 = bl (2.15) 
where bo and cq are constants. Under this normalization, we obtain the equation, 



XYX = — 1 + C ^ Y 

q Y + coqqt ' 



bo n -t 



where the variables X = X{t) and Y = Y(t) are defined by 

X = ~ a -^q^, Y=±q-W , 
bo V c q 



(2.16) 



(2.17) 



respectively, and denotes the value at t — 1. Equation ()2.16j) is essentially the same as 

<?-Pm cm- 

We finally comment on the continuous limit of (/-Pm to Pm- The simplest way is to 
take the limit on the level of equation (|2.14|) : replacing t, K$, Vi and ctj by t/h, 1 + hKi, hVi 
and 1 + (Xih (i = 0, 1), respectively, and taking the limit h — » 0, we obtain equation (|2.5|) . 



2.3 Birational representation of affine Weyl group W"(A^ X 



(1) 



In the following birational representation of (extended) affine Weyl group W(A\ 
A)p\ is presented: 

Theorem 1 (|9j). Let C(a{, fi (i = 0, 1, 2)) be the field of rational functions in a% and f%. 
We define the transformations Sj (i = 0, 1, 2) and ir acting on C(<Zj, fi (i = 0, 1, 2)) by 



Si(aj) = a j a i l \ 
ir(xi) = x i+1 , (x = a,f) 



Si ff ) = f ( ai + h 



(hj = 0, 1, 2) 



i G Z/3Z, 



(2.18) 



respectively, where A = (aij)'fj =0 is the generalized Cartan matrix of type and U 
( u ij)ij=o * s an orientation matrix of the corresponding Dynkin diagram: 

-1 



A 



2 -1 -1 
-1 2 -1 
-1 -1 2 



U 



1 
-1 1 

1 -1 



(2.19) 



Moreover, we define the transformations wq, W\ and r by 

aiO>i+\\ a i-l a i + Oj-l/i + fi-lfi) 



wo(fi) 

Wl(fi) 

r(h) = 



fi-l(aidi+i + aifi + i + fifi+i) 

_ 1 + Ojfi + 

ajOi+l/j+l(l + di-lfi-l + ai-lO>ifi-lfi) ' 

1 

A' 



wo(ai) = toi (ai) = r(di) = a*, (i G Z/3Z). 



(2.20) 
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Then, the transformations (so, s±, s 2 , n, wo, W\, r) generate the extended affine Weyl group 
of typt 
Hons, 



of type A± x A% ■ Namely, these transformations satisfy the following fundamental rela 



s? = l, {s lSi+l f = l, 7T 3 = 1, n Si = S i+1 TT (i = 0,1,2), (2.21) 
Wq = wf = r 2 = 1, rvjQ = W\r, (2.22) 

and the actions of {sq, si, s 2 , tt) and (wo,wi,r) commute with each other. 

One can construct translations Tj (i = 1,2,3,4) by 

Ti = 7rs 2 si, T 2 = siirs 2 , T 3 = s 2 siir, T 4 = rw . (2.23) 

We note that these translations satisfy T{Tj = TfT^ = 1,2,3,4) and T1T2T3 = 1. The 
actions of Tj on and c are given by 



Ti(a ) 


= qao, 


Ti(ai) = g _1 ai, 


Ti(a 2 ) = a 2 , 


Ti(c) = c, 


T 2 (a ) 


= a , 


T 2 (a 1 ) = qa 1 , 


T 2 (a 2 ) = q~ 1 a 2 , 


T 2 (c) = c, 


T 4 (a ) 


= ao, 


T 4 (ai) = ax, 


T 4 (a 2 ) = a 2 , 


T 4 (c) = qc. 



(2.24) 

respectively. We also remark that aoa±a 2 = q and /0/1/2 = °c 2 are invariant with respect 
to the actions of W(A^ x A*ip) and W(A^), respectively. 

The above affine Weyl group actions determine three birational, commutative, dis- 
crete flows on root lattice of type A± x A 2 . Regarding one of the directions as the "time 
evolution", other directions can be considered as those of Backlund (Schlesinger) trans- 
formations. 

In [S], T 4 is chosen as the time evolution, whose explicit action is given as 

™ / , s , 1 + a 2 / 2 + a 2 a f 2 f 1 + a f + a ai/ /i 

Ta{Jo) = aoai/i— 7— 77, r 4 (/i) = a x a 2 f 2 — — — , 

1 + ao/o + aoaiJo/i 1 + 01/1 + 0102/1/2 

T/n 1 + ai/i + a\a 2 f\f 2 , , 

r 4 (/ 2 ) = 0200/0— 5— F-F- ( 2 - 25 ) 

1 + a 2 / 2 + a 2 ao/ 2 /o 

Choosing one of the ^4 2 direction, for example Tj, as the time evolution, we have 
1 + ao/o qc 2 1 + a /o 



Tl(/l)=/ 2 



ao + /o /0/1 ao + /( 







*r'cA)=A£^ = fi£±4. (2.2 6 ) 

1 + ai/i /0/1 1 + ai/i 

Operating T£T%Tg (n, u, N € Z) on equations (l2~23|) and iOHl . denoting T?T%T^(fi) = 
fi(n; v, N) (i = 0, 1, 2), we obtain g-Prv, 

~ „ 1 + a 2 q- p f 2 + a 2 a q n -' / f 2 f 

fo = ao«i<7 h 



fi = aia 2 q n f 2 



1 + a q n f + a aiq u f fx 
1 + aog w /o + aoaiq u fofi 

l + Oi0-"+Vl+OlO 2 g-»/l/2' 



/2 = a 2 a q f — — — , 2.27 

1 + «29 ^/2 + a 2 a g n ^/2/0 

and g-Pm (|l.ll) . respectively. In equation (|2.27|) . denote /j = /i(n; 1/, iV+l) (i = 0, 1, 2). 
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Remark 3. By construction, g-Pm describes one of the Backlund transformations 
of q-Pw (|2.27f) . and vice versa. Moreover, under the specialization ()1.4|) . g-Prv (|2.27|) 
coincides with the Backlund transformation of q-Pn (|1,3[) proposed by Joshi et al jHj. 

g-Pni admits the symmetry described by the affine Weyl group of type x . In 
fact, we can easily verify the following fact: 

Proposition 1. We define w' , w[ by 

w' = sosxso, w[ = s 2 . (2.28) 
Then, w' and w[ satisfy the fundamental relations ofW(A^), 

K) 2 = K) 2 = 1- (2-29) 
Moreover, the actions of (w' , w[) commute with the actions of T\ and (wo,w\,r). 

3 Particular solutions of Riccati type 

3.1 Solutions in terms of Jackson's modified qr-Bessel functions 

Many of Painleve and discrete Painleve equations are known to admit 1-parameter family 
of particular solutions expressible in terms of (g-)classical special functions. Such solutions 
are obtained by finding the special cases where the Painleve or discrete Painleve equations 
are reduced to the (discrete) Riccati equations. 

For example, q-Piy Q2.27JI admits the following Riccati type solutions: 

Proposition 2 (|9j). q-Pjy j2.27[ ) admits particular solutions given by 

sn+MH-l G n(N) . q G n (N + l) 

f °- a ° Cq G n (N + iy h --^T G n (N) ' h ~~ l (3 - 1} 

for v = 0, where G n (N) is a function satisfying the contiguity relations 
G n+1 (N + 1) = G n (N + 1) + a 2 c 2 q 2n+2N+1 G n (N), 

G n+1 (N) = G n (N + 1) + c 2 q 2N+1 G n (N). (3.2) 

In particular, equations (|3.2j) admit polynomial solutions in q N+l l 2 c for ao = q, and the 
polynomials H n (x) = (q N+1 / 2 c)- n G n (N), x = (q N+1 / 2 c + q~ N - 1 / 2 c~ 1 )/2 are identified as 
the continuous q-Hermite polynomials [H]. 

In the above solutions, the variables N and n play a role of independent variable and 
parameter, respectively. Since g-Pni describes a Backlund transformation of q-PrV; the 
above solutions of g-Piv can be also regarded as those for g-Pni only by exchanging the role 
of iV and n. For notational simplicity, we introduce a parameter [i and a function I^in) 
by 



q 



N+1 l 2 c = q», J M (n) = (q N+l ' 2 c)- n G n {N). (3.3) 
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Then equations (|^.2jl are rewritten as 

I^+xin + 1) = <T%+i(n) + agg 2n+ %(n), 

I> + 1) = g~% + i(n) + g%(n), (3.4) 
respectively. We note that the three-term relation of i^(n) with respect to n, 

I> + 2) = + q-") I> + 1) - (l - a 2 q 2n ) I», (3.5) 
is nothing but Jackson's g-modified Bessel equation 
Proposition 3. Let I^{n) be a function satisfying equation {3.$ - Then 

Jo — q ^a - -—, fi — , {A.b) 

safe/y g-P m flip v = and q N+1 / 2 c = q» . 

We remark here that essentially the same solution is also constructed in [3J [2U] . By 
successive application of the Backhand transformation T2 on these solutions, we obtain 
"higher-order" solutions for v G Z>o- Moreover, such solutions admit the following deter- 
minant formula. 

Theorem 2. Let L^iji) be a function satisfying equation {3.$ - For each v € Z>o, we 
define a v x v determinant (fv(n) by 

rf(n) = det (J M (n - i + j))i< M <, 6 Z> ), y#(n) = 1. (3.7) 

Then, 



<pfti(n)(ffi(n) ao $ (n-l)^ +1 (n) 



safe/y q-Pm n~7]J with v € Z> and q N+1 / 2 c = g-^. 

Since Theorem |2] is obtained by the interpretation of the case of g-Piv ;9j , we omit 
the proof. We also note that it is possible to construct similar determinant formula for 
negative v [Oj. 



3.2 q-Pm specific solutions 

In the previous section, we have discussed such solutions that re, v and N (or fj,) played 
a role of independent variable, determinant size and parameter of the solution, respectively. 
It is possible to construct another Riccati type solution in which the roles of v and N (or //) 
are exchanged. 

We first look for the special values of N and c which admit decoupling of g-Pm f|l . ID 
into discrete Riccati equation. In fact, it is easy to verify the following proposition: 
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Proposition 4. q-Pm kl.l\) is decoupled into the discrete Riccati equation, 

fx = f , fi = -aoaiq , (3.9) 

lo/ien iV = and c = 1. 

Eliminating /i from equation (|3.9|) . we obtain a discrete Riccati equation for /o, 

_ „n+l t f„2„2„2i/-l „2_2n+lA j? 

Jo — r~ ^7 • {a. xv) 

1 + a <? n /o 

Standard linearization of equation Q3.1U[) gives the following result: 

Proposition 5. Let H v (n) be a function satisfying the following linear difference equation, 
H u (n + 1) = [1 + ala\q 2u - a 2 q 2n ] H v {n) - a 2 ajq 2u H u (n - 1). (3.11) 

Then, 

1 H v (n + l)-H v (n) 



a q n H v [n) 



/l "" a ° ai9 H v {n + l)-(l-a* q *")H v {n)> (3 ' 12) 

satisfy q-P\u kl.l\) for N = and c = 1. 

One can investigate the evolution with respect to v by substituting the above expres- 
sion into the Backhand transformations of g-Pin described in Section 2.2. After some 
manipulation we arrive at the following result: 

Theorem 3. Let H v {n) be a function satisfying the following contiguity relations, 

H u (n + 1) - H v (n) = -a 2 q 2n H u ^{n), (3.13) 
H u+1 (n) - H u {n) = -a\a 2 q 2v H v (n - 1). (3.14) 

Then, 

n H u _ 1 (n) 1 H v (n) 

/o " ~ aoq ~M^T' h ~ alaxq^H^in-lY ^ 

satisfy q-P\u {HZji with N = and c = 1. 

Remark 4. As is expected from the symmetry described in Section 2.2, the contiguity 
relations (|3.13j) and 1)3.14(1 for H u {n) are symmetric with respect to n and v. In particular, 
the three-term relation with respect to v is obtained as 

ff„+l(n) = [1 + a 2 q 2n - a 2 a 2 q 2u ] H u {n) - a 2 q 2n H v ^(n). (3.16) 
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We obtain higher-order solutions in the form of rational functions in H v (n) by applying 
the Backlund transformation T4 to the above solution. Furthermore, these solutions admit 
the following determinant formula: 

Theorem 4. For each N € !<>o, we define an N x N determinant (f>^ (n) by 

^(n) = det(fT 1 ,+ i _i(n + i-l)) 1 < i)i < JV (N € Z >0 ), 0°(n) = 1. (3.17) 

T/ten, 

. = 1 ^ +1 Nfc(n-l) 

71 a§o ig -+-+^-2 ^+i(n-l)e(") ' 

satisfy q-Pm with c = 1, N £ Z>o- 

We next prove Theorem HJ It is a direct consequence of the following "multiplicative 
formula" with respect to (p. 

Proposition 6. The following formulas holds. 



1 + ao9 n /o = q~ 2N Yv , / ";7v 1 )'\ ~ . ( 3 - 19 ) 



A + 1) 

JV+1 ^n)^_ 1 (n) 



! + Jb_ _ .^^j^MM (3 20) 



l + oi? Ji - _ 2 _,r„_n ,jv+i, 777777 ; > ^.21) 



1 /1 1 e/W^^-l) r , 22 . 

In fact, it is an easy task to derive g-Pm by using the above multiplicative formulas 
and equation ()3.18j) . Furthermore, Proposition follows from the following proposition. 

Proposition 7. <ft^ (n) satisfies the following bilinear difference equations. 

^ +1 (n)^_x(n) - a 2 oq ^ +N U^(n)^(n) 

= q- 2N ^ +1 (n+l)^_ 1 {n-l), (3.23) 
q 2N ^+\n)^(n) - aga?^"" 1 )^* - l)^(n + 1) 

= ^ +1 (n)$Li(n), (3.24) 
og^^C" " l)^(n) - T™tf! + \n)4$-x{n - 1) 

= ^ +1 (n-l)^_ 1 (n), (3.25) 

d 2{v - l) m{n ~ l)^(n) - g-^+V^lC" " 1) 

= ^ r i 1 (n)^ r (n-l). (3.26) 

For example, equation ()3.19j) is derived from equation ()3.23(l by dividing both sides by 
4>^ +l (n)4>^_ 1 (n). Therefore, the proof of Theorem 0] is now reduced to that of Proposi- 
tion [7| We will give it in the next section. 
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3.3 Proof of Proposition 

Our basic idea for proving Proposition is the same as that was done in [§] . Bilinear 
difference equations are derived from the Plucker relations, which are quadratic identities 
among determinants whose columns are shifted. Therefore, we first construct such "dif- 
ference formulas" that relate "shifted determinants" and (j)^ (n), by using the contiguity 
relations of H u (n). We then derive bilinear difference equations with the aid of difference 
formulas from proper Plucker relations. We take equation ([3.24)1 as an example to show 
this procedure explicitly. For other equations, see Appendix A. 

We introduce a notation, 



n 



H u (n) 
HJn + 1) 



H v+ i{n) 
H u+1 (n + l) 



H u {n + N-l) H u+1 (n + N-l) 
.JV-1J, 



Hv+N-i(n) 
H v+N -\{n + 1) 

H v+N - 1 (n + N- 1) 



(3.27) 



where k n denotes a column vector, 



I H u+k (n) 
H v+k {n+l) 



(3.28) 



We note that we abbreviate the suffix n when unnecessary. 
We next construct the difference formula. 

Lemma 1 (Difference formula I). 



|0,1, ■•■ ,JV-l| = ^(n), 

|0 n+1 , 0, • • • , N - 2| = (-alaiy {N - 2) g -(^-a>(*-i)^ (n + 1}j 



|ln+l, 0, 



,N 



(-Oo°i) 



2 2\~(N-2) _(2 l /+JV-2)(JV-l) iN 



-(n + 1). 



(3.29) 
(3.30) 
(3.31) 



Proof. Equation (|3.29|) is nothing but the definition equation ([3.27)1 . Subtracting the 
(k — l)-th column from the k-th column of the left hand side of equation ([3.29)1 and using 
the contiguity relation ([3.14)1 for k = N, N — 1, . . . , 2, we have 



n 



H u+N ^ 2 (n) 
H v+N ^ 2 (n + 1) 



H v+N -x{n) - H u+N - 2 (n) 
H u+N ^i{n + 1) - H u+N ^ 2 (n + 1) 



H v+N „ 2 (n + N-1) H v+N ^{n + N — 1) — H u+N „ 2 (n + N - 1) 
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H u+N - 2 {n) 
H v+N _ 2 (n + 1) 



-ala\q^+ N -VH u+N _ 2 {n-l) 
-ala\ q ^+ N -^H v+N _ 2 {n) 



H u+N - 2 (n + N-1) -a 2 a 2 iq 2 ^ +N ~^H u+N ^ 2 (n + N - 2) 



-ala\q^ +N -V 

H v (n) 
H v {n + l) 



H v+N - 2 (n) 
H u+N _ 2 (n + 1) 



H v+N ^ 2 {n - 1) 
H u+N ^ 2 (n) 



H v {n + N-l) ••• H u+N ^ 2 (n + N - 1) H v+N _ 2 (n + TV - 2) 



2 2^-1 „(2iH-JV-2)(iV-l) 



H v {n) 
H v {n + 1) 



H v (n - 1) 
HJn) 



H u+N _ 2 (n - 1) 
H u+N ^ 2 (n) 



H u {n + N-l) H u (n + N-2) ••• ^^(n + iV - 2) 



(3.32) 



which is nothing but equation (|3,3U|) , Moreover, subtracting the second column multiplied 
by a 2 ) a\q 2u from the first column of the right hand side of equation ()3.32j) and using the 
contiguity relation (|3.14|) . we obtain 



n 



H u+1 (n) H v (n-1) 
H v+1 (n+l) H v {n) 



H v+1 {n + N-l) H v {n + N-2) ■ 
which is nothing but equation (|3.31|) . 
Now consider the Pliicker relation, 

= |l„+i,l,--- ,N-1,N\ x |0, 1, - - - ,N-l,(p 2 \ 

- |0, 1, ,JV-1,JV| x |l n+1 ,l,l,--- , N — l,ip 2 \ 

- ,N-l,N,(p2\ x |l n+ i,0,l,-- - ,JV-1|, 

where <p 2 is the column vector, 



H u+N - 2 (n - 1) 
H v+N _ 2 (n) 

H v+N - 2 {n + N - 2) 



(3.33) 



(3.34) 



<f2 



( 1 \ 



v o y 



(3.35) 



We obtain equation (|3.24j) by applying Lemma ^ to equation ()3.34j) after expanding with 
respect to the column <p 2 . This completes the proof of bilinear equation (|3.24[) . 
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A Difference formulas and Pliicker relations 



In this appendix, we provide with data which are necessary for the proof of Proposition [Tj 
We first give the difference formulas. 

Lemma 2 (Difference formula II). 

|0', 1', • - • ,N-l'\ = <tf?{n), (A.l) 

10^, 0V-. ,N-2'\ = (-aiy {N - 2) q-^-^U^An), (A.2) 

,N-2'\ = (-aiy {N - 2) q-^-W-VtfL^n), (A.3) 



r rv 



where 



( H u {n + k) 
H u+ i(n + k) 



(A.4) 



V 



We next rewrite <ft^ (n) in terms of Toeplitz determinant with respect to n by using the 



contiguity relations Q3.13JI and (|3. 14f) as 

H v (n+l) 



2„2\ N ( N - 1 )/ 2 N(N-l)(3v+N-2)/3 



H v (n) 
H v {n-l) 



HJn) 



H v {n + N-l) 
H u (n + N-2) 



H v (n-N + l) H u (n-N + 2) ■■■ H u (n) 

= (- a 2 a 2)^-l)/2 9 iV(iV-l)(3,+iV-2)/3 (n) 

Then we obtain the third difference formula: 
Lemma 3 (Difference formula III). 



|0,1, 

|^ + i,o,- •• ,W 



\l v+ i,0, ••• ,N 



where 



k, 



( H u (n + k) \ 
H v (n + k-l) 
H u (n + k-2) 



\ 



( H u (n + k) \ 
q- 2 H u (n + k-l) 
q- 4 H u (n + k-2) 



J 



\ 



(A.5) 



(A.6) 
(A.7) 
(A.8) 



(A.9) 



/ 



Lemmas 13 and El are proved by applying the same procedure as that for Lemmadon 
the determinants, 



H v {n) 
H u+ i(n) 



H v {n + 1) 
H v+l {n + l) 



H v+N _i(n) H v+N -i{n + 1) 



H v {n + N-l) 
H v+l {n + N - 1) 

H u+N ^{n + N -1) 



(A.10) 
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and 



H u (n) 
H v (n-1) 



H u {n) 



H u (n-N + 1) H v {n-N + 2) 



H u (n + N-1) 
H u (n + N-2) 

HJn) 



(A.ll) 



respectively, by using the contiguity relation ()3.13j) . 

The bilinear equations (|3,23|) - ()3.26|) are derived by considering proper Plucker relations. 
We finally give the list the Plucker relations and difference formulas which are necessary 
for the derivations. This completes the proof of Proposition [7| and thus Theorem |1] 

Equation 
Plucker relation 



0= Ih-i,!,- •• ,N-1,N\ x 0,1,-.. 



|0,1,-.. ,N-1,N\ x IIh.1,1, -•• ,N-l,<p 2 \ 



- ,N-l,N,ip 2 \ x 0,1,- ■■ ,N-1\. 

Difference formula Lemma EH 

Equation (|3.25|) : 
Plucker relation 

0= .JV-l'.JV'l x 10', 1', - - - .JV-1',^1 

- |0', 1', • • ■ ,N-1',N'\ x |1'„ +1 ,1',1',.-- 
-|1',-.. ,JV-l',JV',^i| x ll^O'.l',.-- ,JV-1'|. 

Difference formula Lemma [51 

Equation <$57Z6§ : 
Plucker relation 

0= |l n+ i,l,--. ,JV-1,JV| x |0,1,--. ,N-l,(p!\ 

- |0, 1, ,N- 1,N\ x |l n+1 ,l,l,-.. ,N — l,<pi\ 

- |1,-.. ,N-l,N,ipx\ x |l n+1 ,0,l,--- ,N-1\. 

Difference formula Lemma ^ 

In the above data, ipi and </?2 are column vectors given by 





V i / 



¥2 



( 1 \ 





(A.12) 
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